Stochastic resonance in the driven Ising model on small-world networks 
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We investigate the stochastic resonance phenomena in the field-driven Ising model on small-world 
networks. The response of the magnetization to an oscillating magnetic field is examined by means 
of Monte Carlo dynamic simulations, with the rewiring probability varied. At any finite value of 
the rewiring probability, the system is found to undergo a dynamic phase transition at a finite 
temperature, giving rise to double resonance peaks. While the peak in the ferromagnetic phase 
grows with the rewiring probability, that in the paramagnetic phase tends to reduce, indicating 
opposite effects of the long-range interactions on the resonance in the two phases. 

PACS numbers: 89.75.Hc, 05.40.-a, 75.10.Hk 
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I. INTRODUCTION 



It has been known that a periodically modulated 
bistable system may display stochastic resonance (SR), 
arising from the cooperative interaction between random 
noise and periodic modulation Namely, when a sys- 
tem with an energetic activation barrier is subject to pe- 
riodic but weak driving, the inherent thermal stochastic 
noise can enhance the signal out of the system rather 
than weaken it. Such SR phenomena, which have various 
practical applications have been extensively investi- 
gated in a variety of systems: In parallel to studies of 
systems with relatively small numbers of degrees of free- 
dom, the research focus has been shifting to extended 
complex systems with many degrees of freedom, where 
interesting collective dynamic behavior can emerge ||, [| . 
In those studies of SR in extended systems, the underly- 
ing connection topology of dynamic variables has usually 
been assumed to be regular. However, recent studies of 
computer networks, neuronal networks, biochemical net- 
works, and even social networks, have revealed that many 
real systems in nature possess quite complex structures, 
which can be described neither by regular networks nor 
by completely random networks || . 

In this paper we consider an Ising model on Watts and 
Strogatz (WS) type small-world networks [[| and study 
the SR behavior in the presence of temporally oscillating 
external magnetic fields. The WS network is character- 
ized by a short characteristic path length and high clus- 
tering j(| , both of which are commonly observed proper- 
ties of real networks in nature. Accordingly, we believe 
that the study of SR in an extended system can be made 
more realistic if one uses the WS network as the under- 
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lying topology. 

There are five sections in this paper: Section II intro- 
duces the system, i.e., the Ising model, driven by oscil- 
lating magnetic fields, on small-world networks. Defined 
in Sec. Ill is the dynamic order parameter, which de- 
scribes conveniently the dynamic phase transition. The 
relaxation behavior is examined and the relaxation time 
is computed. Section IV is devoted to the stochastic res- 
onance phenomena, characterized by double resonance 
peaks in the occupation ratio. The effects of long-range 
interactions on the resonance are investigated, revealing 
opposite trends in the ferromagnetic and paramagnetic 
phases. Finally, the main result is summarized and dis- 
cussed in Sec. V. 



II. ISING MODEL ON SMALL- WORLD 
NETWORKS 



The WS network in the present paper is constructed 
following Ref. Q: First, a one-dimensional regular net- 
work with only local connections (of range k) is con- 
structed with the periodic boundary condition. Next, 
each local link is visited once, and with the rewiring prob- 
ability P removed and reconnected to a randomly chosen 
node. After the whole sweep of the entire network, the 
average number of shortcuts in the network of size N is 
given by NPk. Throughout this paper, the interaction 
range k is set equal to two for convenience; longer ranges 
(k > 2) are not expected to lead to any qualitative dif- 
ference. After the WS network is built as above, an Ising 
spin is put on every node, and an edge (or a link) con- 
necting two nodes is regarded as the coupling between 
the two spins at the two nodes. 

The Hamiltonian for the field-driven Ising model on 
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the WS network thus reads 



(i) 



where J is the coupling strength, ct; (= ±1) is the Ising 
spin at node i, the neighborhood Aj of i stands for the set 
of nodes connected to i (via either local edges or short- 
cuts), and an oscillating magnetic field h(t) = ho cos fit 
is applied with the driving amplitude ho and frequency 
O. We perform Monte Carlo (MC) dynamic simulations, 
employing the heat bath algorithm J7j and measuring the 
time t in units of the MC time step. For thermalization, 
we start from sufficiently high temperatures and lower 
the temperature T slowly with the increment AT = 0.02 
(in units of J/ks with the Boltzmann constant fee). The 
driving amplitude and frequency are mostly taken to be 
ho = 0.1 and £1 = 0.1 although different frequencies are 
also considered. While simulations are performed at a 
given temperature, the data from the first 4 x 10 4 MC 
steps are discarded, which turns out to be sufficient for 
stationarity, and measurements are made for next 10 4 
MC steps. Networks of various sizes, up to N = 6400, 
are constructed as described above, and averages are per- 
formed over 100 different network realizations. 

In the absence of the long-range interaction (P = 
0), the network structure reduces to that of the one- 
dimensional regular network with only local couplings. 
Accordingly, when P = 0, the driven Ising model de- 
scribed by Eq. (|l|) as well as the undriven model (ho = 0) 
should not exhibit long-range order at finite tempera- 
tures. For P ^ 0, on the other hand, it has been found 
that the (undriven) Ising model displays ferromagnetic 
order at finite temperatures |8). This suggests that the 
driven Ising model in Eq. (|l|) should undergo a dynamic 
phase transition at a finite temperature unless h(t) is too 
large. When all nodes are fully connected, Eq. ([!]) de- 
scribes the infinite-range Ising model, where double SR 
peaks have been observed and argued to be a generic 
property of the system with a continuous dynamic phase 
transition || . In view of this, similar double SR peaks are 
naturally expected when P / 0; this contrasts with the 
case P = 0, where the absence of a dynamic phase tran- 
sition implies the emergence of just a single SR peak p0|. 



III. DYNAMIC PHASE TRANSITION 

The dynamic phase transition in the system can be 
conveniently described by the dynamic order parameter. 
We first measure the magnetization 
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FIG. 1: Dynamic order parameter Q versus the temperature 
T in the WS network of size N = 200, 400, 800, 1600, 3200, 
and 6400, at the rewiring probability P = 0.5. Inset: Dy- 
namic fourth-order cumulant Un has a unique crossing point 
at T c » 3.13 (in units of J/ks)- 



and take the time average of m(t) during the n-th period 
of h(t), to obtain 



O 
2^ 



t„+i 



dt m(t) 



(3) 



tn 



with t n = 2irn/Q. The dynamic order parameter is then 
defined to be 



Q = lim Q n , 



(4) 



which takes different values below and above the criti- 
cal temperature T c : Figure [j] shows that Q ^ in the 
ferromagnetic phase at T < T c and Q — in the para- 
magnetic phase at T > T c , with T c w 3 for P = 0.5. 
In order to determine T c more precisely, we measure the 
dynamic fourth-order cumulant [O, rL2| 
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where (• • • ) and [• • •] denote the time average and the 
average over different network realizations, respectively. 
At T = T c , the cumulant Un should have a unique value 
regardless of the size of the system; this yields the esti- 
mation T c w 3.13, as shown in the inset of Fig. 0. 

To investigate the time relaxation behavior, we be- 
gin with the initial condition m(t=0) = 1 (i.e., a- L = 1 
for all i), and measure Q n in Eq. (pi) as a function of 
t n during MC simulations. As time proceeds, Q n ap- 
proaches the dynamic order parameter [see Eq. (Q)] with 
the value Q = and Q ^ in the high- and low- 
temperature phases, respectively. We find that the re- 
laxation of Q n is very well described by the exponential 
form, Qn — Q oc e~*"/ T with t n = 27rn/f2, which defines 
the relaxation time r |l]| . Figure || shows the relaxation 
behavior for P — 0.5 and 0.8 at T = 3.5; it is observed 
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FIG. 2: Relaxation of the dynamic order parameter. Q n — Q 
versus t n at temperature T = 3.5 is plotted in the semilog 
scale for the rewiring probability P — 0.5 and 0.8. As the 
measurement time t n is increased, Q n is shown to approach its 
steady state value Q. Dashed lines represent the least-square 
fit to the exponential-decay form, Q„ — Q cx e _t "^ r , from 
which one can obtain the relaxation time r. Inset: Inverse of 
the relaxation time 1/r versus the temperature T. 



that the exponential decay form describes the simula- 
tion results very well. We repeat the same procedure at 
other temperatures and obtain the temperature depen- 
dence of the relaxation time, shown in the inset of Fig. ||. 
As expected from the existence of the finite-temperature 
dynamic phase transition, the relaxation time diverges 
near the dynamic transition point (T c w 3.13 and 3.20 
for P — 0.5 and 0.8, respectively), and decreases as we 
move away from T c in both directions || . 



IV. STOCHASTIC RESONANCE 

In this section we study how the long-range interac- 
tions influence SR by varying the rewiring probability P. 
The SR behavior is conveniently captured by the occu- 
pancy ratio R, first introduced in Ref. || and defined to 
be the average fraction of the spins in the direction of the 
external field Q : 

j£ — / num ber of spins in the direction of h(t) \ 
\ total number of spins / 

In other words, R measures how many spins follow the 
oscillating magnetic field. It is easy to understood that 
R has the value 1/2 in both low- and high-temperature 
limits (see, e.g., Ref. ||) and becomes larger near the 
SR temperature, reflecting that more spins follow the 
external driving. 

In Fig. ||, R in the system with P = 0.5 is plotted as a 
function of the temperature for various sizes N. The dou- 
ble SR peaks are clearly exhibited, with the dip emerging 
in the vicinity of the dynamic phase transition temper- 
ature T c determined from the crossing of U n . As the 



FIG. 3: Size effects of the occupancy ratio R versus the tem- 
perature T for P — 0.5. The finite-size effects are stronger for 
the SR peak at the lower resonance temperature but become 
negligible for N > 800. Insets: Size-dependence of the SR 
temperature T§x and of the height of the first resonance 
peak. 

size N is increased, the first SR peak at the lower reso- 
nance temperature Tgl becomes sharper and apparently 
saturates for N > 800, displaying negligible finite-size 
effects. The SR temperature Tg^ and the height 
of the first resonance peak are displayed in the insets 
of Fig. ||, manifesting that there are no substantial size 
effects for N > 800. For the second SR peak at the 

(2) 

higher resonance temperature Tg^, finite-size effects are 

(2) 

insignificant in both the height hg^ and the SR temper- 

(2) 

ature T s ^. We thus conclude that the double SR peaks 
are not merely finite-size effects, and obtain T^] w 2.70 

(2) 

and Tgji w 4.46 from the positions of the two SR peaks 
in Fig. 

In Ref. ||, the positions of SR peaks in the infinite- 
range Ising model have been analytically obtained from 
the time-scaling matching condition: The intrinsic time 
scale given by the relaxation time r should match the ex- 
trinsic time scale of the external driving. Furthermore, 
since r diverges only at T c and decreases as the tempera- 
ture T is raised or lowered from T Cl the presence of double 
SR peaks should be a general property of the system with 
a finite-temperature continuous phase transition. How- 
ever, the specific form of the matching condition for the 
infinite-range Ising model, t = 1 + vT+ ^ _2 , found in 
Ref. H may not hold for the Ising model on the WS 
network studied in this work. Here we reveal the rela- 
tion between the positions of double SR peaks and the 
relaxation time in the following manner: We first draw 

a vertical line at T = (w 4.46) in the inset of Fig. § 
and locate the crossing point with the plot for 1/r, from 
which a horizontal line is drawn. This horizontal line 
corresponds to the time-scale matching condition and its 
crossing point with the lower temperature branch of the 
1 /r plot then yields the position of the first SR peak Tg 1 ^. 
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FIG. 4: (a) Occupancy ratio R versus the temperature T and 
the rewiring probability P in the system of size N = 800. 

(b) Lower and upper resonance temperatures Tg 1 ^ and T s 2 ^ 
as well as the dynamic transition temperature T c versus P. 

(c) Heights hg 1 ^ and of the resonance peaks versus P. 

This gives w 2.70, which is in excellent agreement 
with the value obtained from the first peak of R in Fig. [|. 

The behavior of i? with the temperature T at var- 
ious values of the rewiring probability P is shown in 
Fig. ||. When long-range interactions are not present 
(P = 0), the Ising model does not establish ferromagnetic 
order at finite temperatures, thus not displaying a finite- 
temperature phase transition. Then expected is a single 



SR peak instead of two peaks |9| , which is confirmed in 
Fig. |J. As P is increased, the system behaves more sim- 
ilarly to the mean-field system, and exhibits clear dou- 
ble SR peaks, which, in the mean-field (infinite-range) 
limit, have equal heights Note that the maximally 
connected case (P = 1) of the WS network has O(N) 
long-range connections while the infinite-range system 
has 0(N 2 ) connections. This explains why the double 
peaks of R at P — 1 do not have equal heights in Fig. |^. 

The SR peaks in Fig. || may be located from the con- 
dition dR/dT — except for the first resonance peak at 
P = 0.1, where we have used the condition that dR/dT 
has a minimum positive value. In the insets of Fig. |], the 
obtained SR temperatures and peak heights are shown as 
functions of P. It is noteworthy that as P is increased 
the height of the first peak also grows, whereas hc^ 
reduces. In other words, addition of long-range interac- 
tions tends to enhance SR in the ferromagnetic phase, 
while the same addition suppresses SR in the param- 
agnetic phase. To understand these conflicting effects, 
which appear somewhat counterintuitive, we notice that 
as P is increased, the critical temperature T c shifts to- 
ward higher temperatures, reducing the relative (tem- 
perature) distance to Ta^ in comparison with that to 
Tgl. In the paramagnetic phase, the correlation length 

(2) 

should increase with P since becomes closer to T c ; 
this leads more spins to be correlated and thus resistive 
to the change of the external driving field, resulting in 
the suppression of SR. In the ferromagnetic phase, on 
the other hand, a larger number of shortcuts (i.e., long- 
range connections) increases the distance to T c , yielding 
a shorter correlation length. Accordingly, the spins are 
allowed to follow better the external driving, and the 
SR is enhanced. It is also to be noted that SR in the 
present system tends to saturate for P > 0.7 and that 
almost the same SR as the random network (P = 1.0) 
can be achieved with relatively small numbers of short- 
cuts. Similar enhancing and saturating tendency has also 
been found in the synchronization behavior of the cou- 
pled oscillators on small-world networks j^] as well as in 
the standard (single-peak) SR behavior of the system of 
coupled bistable elements Q . 

We have also considered different driving frequencies 
and investigated how the occupancy ratio behaves. Fig- 
ure H displays the behavior of R depending on the tem- 
perature T for various driving frequencies. The inter- 
val between the two resonance peaks is observed to grow 
as the driving frequency is increased, which is consis- 
tent with the obtained relation between the relaxation 
time and the driving frequency. On the other hand, the 
heights of the two peaks tend to reduce with the driv- 
ing frequency; eventually, the lower peak in the ferro- 
magnetic phase disappears, which reflects that at high 
driving frequencies there exists only one crossing point 
between the relaxation time and the driving frequency. 
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FIG. 5: Occupancy ratio R versus the temperature T for the 
rewiring probability P — 0.3 and several driving frequencies. 
While the interval between the two peaks tends to grow with 
the driving frequency £1, the heights of the peaks in general 
reduce with Q. At high frequencies (Q = 0.3), the lower peak 
in the ferromagnetic phase is observed to disappear. 

V. SUMMARY 

We have examined the stochastic resonance phenom- 
ena in the Ising model, driven by oscillating magnetic 
fields, on small-world networks. Double resonance peaks 
have been found to develop as the rewiring probabil- 
ity is increased, and explained in terms of the diverging 
time scale at the dynamic phase transition. It has been 
demonstrated that the resonance behavior essentially the 
same as that of the random network can be achieved 
with relatively small numbers of shortcuts. The most 



interesting finding in our system is that while stochas- 
tic resonance is enhanced with the number of long-range 
connections in the ferromagnetic phase, it is suppressed 
in the paramagnetic phase. As an interpretation based 
on analogy, we suggest to consider the opinion forma- 
tion in a social system: an instructor and a group of 
students. More long-range connections in this analogy 
correspond to stronger interactions, i.e., more active dis- 
cussions, between students far away, and thus help the 
group behave as a whole with a majority of students hav- 
ing the same opinion (this situation is closely related to 
the enhanced synchronization for larger P in Ref. jL5|). 
However, this opinion developed from active discussions 
among students may not necessarily be the opinion pre- 
ferred by the instructor of the class. The enhanced SR 
peak in the ferromagnetic phase corresponds to the sit- 
uation that more interactions among students help the 
class follow the instructor, with the students having al- 
ready an identical opinion (whatever it is). On the other 
hand, the suppression of SR in the paramagnetic phase is 
interpreted as the situation that more interactions make 
the students, among whom an agreement is not reached, 
not so obedient to the instructor. 
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